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5?R0CEJ)B DE SffiCOBTSTRtTCTION D'DNE ISC^GE X P&SKFJR D'OSI 5HU 
1XE PROJECTIIONS PAR APPLICATION D'UNE TR&ETSFORMEE SN * ' 

ONDSI^ETTE 

DESCRIPTION 

Le sujet de cette invention est un pracSdS 
reconstruction d'une image £ parti r de jeux de 
project ions de cette image et en appliquant une 
transformation en ondelette- 

Les precedes de tomographic consistent a 
examiner un objet inanimfi ou un etre vivant par un 
rSaeau de dStectasUrs mobiles qui prend une aSrie -de 
vues en toumant aufeour de lux. Ces vues sent dea 
projections de la propri&te permettant d'escpriTue*- 
1 1 image (normalement: une absorption d ' un rayonneiaent 
travearsant !• objet on une scintillation d'un cor^a 
radioactif ing£r£ par l 1 objet), e'est-i-dire das sotnmes 
de la propriStS le long de lignes traversant I 1 objet jpt 
dgfinies par la CQllimation des dStecteurs* Chagiie 
d^tecteur manure une projection de 1 1 image k chaqSe 
vue. Quand un - nowibre suffisant de vues et de 
projections a &z& pris, on procSde a l f inversion des 
rSsultats pour obtenir' la valeur de la propriStS a 
cbaque point de 1 'objet ; CQtte inversion est 
comparable Sl l f inversion d<un syst^me d T Equations de 
grandes dimensions et peut s'effectuer ouverteraent par 
dea m$thodas algSbriques, ou, plus f rSquemment r par dec 
tnSthodes analytiqueq par lesquelles on applique des 
operations numeriques successives auat projections, ssns 
inverser direefcsmant le systdme. Un grand nombre de 
" mSthodes existent, parmi lesq^ielle? on citejra celle qui 
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PROCBDE DE RECONSTRUCTION D'UNE IMAGE A PARTIR D'UN JEU 
DE PROJECTIONS PAR APPLICATION D'UNE TRANSFORMEE EN 

ONDELETTE 



DESCRIPTION 

Le sujet de cette invention est un procede 
de reconstruction d'une image a partir de jeux de 
projections de cette image et en appliquant une 
transformation en ondelette. 

Les procedes de tomographie consistent a 
examiner un objet inanime ou un etre vivant par un 
reseau de detecteurs mobiles qui prend une serie de 
yues en tournant autour de lui . Ces vues sont des 
projections de la propriete permettant d'exprimer 
1' image (normalement une absorption d'un rayonnement 
traversant l'objet ou une scintillation d'un corps 
radioactif ingere par l'objet), c'est-a-dire des sommes 
de la propriete le long de lignes traversant l'objet et 
definies par la collimation des detecteurs. Chaque 
detecteur mesure une projection de 1 ' image a chaque 
vue. Quand un nombre suffisant de vues et de 
projections a ete pris, on procede a 1" inversion des 
result at s pour obtenir la valeur de la propriete a 
chaque point de l'objet ; cette inversion est 
comparable a 1 • inversion d'un systeme d' equations de 
grandes dimensions et peut s'effectuer ouvertement par 
des methodes algebriques, ou, plus frequemment, par des 
methodes analytiques par lesquelles on applique des 
operations numeriques successives aux projections sans 
inverser directement le systeme. Un grand nombre de 
methodes existent, parmi lesquelles on citera celle qui 
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est d£crite dans le brevet fran9ais 2 615 619 qui est 
le premier de cette £quipe de recherche, et le brevet 
£ran<?ais 2 810 141 plus recent qui present e quelquee 
re0^etQbl«mces avec le proc6d§ qui sera d£crit icx* On 
pent aujaei citer un article {<3©2»ce» — JcS — g&ference— 
ceng)3^te) — r^de.. grass. > Ces operations aboutiasent A ce 
^ qu T on appelle la rg^roproj action dee mesnres, c'est-a- 
^ dire an calcul de la valeur de la propriety pris^ & 

if lp des points de \a ligae de la projection. 

^ "70#fc 10 11 P eut ^tre avantageux de travailler avec 

des rSsultat© e^primea dans le domaine frequent i el de 
Fourier r comma le prouve le deroci&na document. Des 
transformers numeriques d 1 autre nature ont aussi St£ 
employees. 

IS On reportera aux figures 6, 7 et S poire* 

•una description concrete et sclx&fnatiqua d'un proc€d£ de 
prise de mesqres • tomographiques . Utae source P de 
rayonnement et un syst&me de detection 3 sont mobiles 
le long d<un cadre* axmulaire 2 a des positions 
20 dia^tralement qpposSes, et mx rayonnement conique 
originaire de la source P atteint Is syst&ae de 
detection 3 apr&s avoir traverse I'dbjet 1 k 6tudier« 
lia piSce essentielle du • systSme de detection 3 est un 
r€seau 4 bidimensiomsel de d&tecteurs 5. Des 
2S projections R de 1 1 image tridiraensionnelle de I'bbjet 1 
sont mesurees par ceux dee d£tecteurs S qui sont inclns 
dans un potprtour 15 de 1' "ombre" de I'objet 1* XJn 
notnbre important de vuee de ce genre est pris 3, autant 
d 1 angles 0 diff brents <i 11 orientation du reseau 4 des 
. 30 dStec tenors S. 11 est frSquemrnent considSre, 
convent ionnel lenient pour les calcule, un r€seaxx fictif 
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est d6crite dans le brevet fran<?ais 2 615 619 qui est 
le premier de cette 6quipe de recherche, et le brevet 
francpais 2 810 141 plus recent qui pr6sente quelques 
ressemblances avec le proc€d6 qui sera d€crit ici. On 
peut aussi citer un article (donnez la reference 
complete) ... de Grass. Ces operations aboutissent a ce 
qu'on appelle la r^troprojection des mesures, c'est-a- 
dire au calcul de la valeur de la propriety prise a 
chacun des points de la ligne de la projection. 

II peut etre avantageux de travailler avec 
des resultats exprimes dans le domaine frequentiel de 
Fourier, comme le prouve le deuxieme document. Des 
transformees num£riques d f autre nature ont aussi ete 
.employees. 

On se reportera aux figures 6, 7 et 8 pour 
une description concrete et schematique d'un proc^de de 
prise de mesures tomographiques . Une source F de 
rayonnement et un systeme de detection 3 sont mobiles 
le long d'un cadre annulaire 2 a des positions 
diametralement opposees, et un rayonnement conique 
originaire de la source F atteint le systeme de 
detection 3 apres avoir traverse I'objet 1 a gtudier. 
La piece essentielle du systeme de detection 3 est un 
rlseau 4 bidimensionnel de d£tecteurs 5. Des 
projections R de 1 1 image tridimensional le de I'objet 1 
sont mesurees.par ceux des detecteurs 5 qui sont inclus 
dans un pourtour 15 de 1 ■ " ombre " de 1 1 ob j et 1 . Un 
nombre important de vues de ce genre est pris a autant 
d' angles 9 differents d ! orientation du reseau 4 des 
d gtecteurs 5. II est frgquemment considers, 
convent ionnellement pour les calculs, un reseau fictif 
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4 f de d&tecteur 5* sur un plan de detection 5det qui 
est parall&le au r£seau 4 rfiel et passe par le centre O 

du cadre .2. Des coordonnSes p et q sont d€f iniee pour ^> 

* — 

identifier les dStecteurs ^^^d'apres leurs lignes et 
5 leurs eolonnes. Des calculs de rearrangement , courants 
dans I 1 art, permetxent de transposer lea assures d'un 
r£seau 4 quelconque au reseau fictifi 4' et d'appliquer 
les aigorithraes de reconstruction a celui-ci f 

brevet frangais 2 615 619 sera rappel€ 
10 ici puisqu'il expliquait en dStail quelles operations 
numeriquep, notairanent celles qui sont dites de filtrage 
et de r$troprojectibn/ permettent d'.obtenir l 1 image de 
1'objet a partir de* ses projections ; des mStfaodes 
diffSrentes existent ftependant. 

-\ 

15 I»e problSme auquel on s B attaque ici est la 

reduction du teiups ou du volume des calculs pour r 

inverser le syst&me des ciesures. 

On saifc qu'il s'agit d'une des limitations 

les plus graves dans les ps:oc£d£s de tomograpbie, et de 
20 nombreuae prooSdSs nouveaux ant StS congua pour le tuSrae 

object if que l f invention, dont celui du deuxiSme brevet 

cit£- 

L*id€e exposSe ici est d' exploiter les 
propriStfis particull^res d'une transformee numSrique, 
25 dite decomposition eri ondel^J^i, des , P r ^^^^ ns P olir r 



pzfivoir des emplacements de pr^^rrax^. nSgligeables ou 
insignifiantes et de ne pas appliquer les ealculs 
d' inversion & ces emplacements , Des perfectionnemenfcs 
off rent encore une plus grande reduction des calculs, 
30 Sous sa forme la plus g&iSrale, 1* invention 

edncexxte un pxccSdS de reconstruction d*une image 
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4 1 de detect eurs 5 ' sur un plan de detection Pdet qui 
est parallele au reseau 4 reel et passe par le centre 0 
du cadre 2. Des coordonnees p et q sont definies pour 
identifier les detecteurs 13 d'apres leurs lignes et 
leurs colonnes. Des calculs de rearrangement, courants 
dans l'art, permettent de transposer les mesures d'un 
reseau 4 quelconque au reseau fictif 4' et d'appliquer 
les algorithmes de reconstruction a celui-ci. 

Le brevet francais 2 615 619 sera rappele 
ici puisqu'il expliquait en detail quelles operations 
numSriques, notamment celles qui sont dites de filtrage 
et de retroprojection, permettent d'obtenir 1 ' image de 
l'objet a partir de ses projections ,- des mSthodes 
diffSrentes existent cependant. 

Le problSme auquel on s'attaque ici est la 
reduction du temps ou du volume des calculs pour 
inverser le systSme des mesures. 

On sait qu'il s'agit d'une des limitations 
les plus graves dans les proced€s de tomographie, et, de 
nombreux procedes nouveaux ont ete con?us pour le meme 
object if que 1' invention, dont celui du deuxieme brevet 
citS . 

L'idee exposee ici est d 1 exploiter les 
propri6tes particulieres d'une transformee numerique, 
dite decomposition en ondelette, des projections pour 
prevoir des emplacements de projections negligeables ou 
insignif iantes et de ne pas appliquer les calculs 
d' inversion k ces emplacements. Des perf ectionnements 
off rent encore une plus grande reduction des calculs. 

Sous sa forme la plus genSrale, 1' invention 
concerne un procede de reconstruction d'une image 
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A ^^^^ 

comprenant na sSrie de d£coroposit:xaiia(^uc^ 

fci or^ft ea , ondelette dormant {Qtesafaa gefe t e ie 

I w ,, i. n M . w jN sy f.- - > 

der pjoJeeLlau ^J^iup^ approximation et 

des siries ^raccessives d$ details nom^^aes de ch&que 
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a-0 



15 



20 



25 



30 




fcgoprgjggj^ jeux de 



jeu, puis de 

projectiox^^fiT imagettes de 1" image £ reconstruire et 
une combinaisoxx dee' iiaagettes a re const mi re par une 
inversion de decomposition en ondelette dormant ladite 
image, caract£ris£ en ce qu'il eomprend, dans chacune 
d«© series et successivement pour des iinagefctes de 
details ayant un noiribre croissant de points/ une 
recherche de portions isusignifiantes, -JugSes pour 
depourvuesr de contezm, et Tine recherche de portions 
homoth&tiques des portions ins ignifi antes pour chacune 
des iinagettes de details qui suivent dams la sfirie ; et 
en ce qu 1 il^jjo^rend -uiier-^mi^&dUDn^des rfitroprojectiox^s 



pour tontee^ K ^ po rti ono ingign jb ficmL^ et CDnfcHgp4bgs^ ^ 
f-pos rtiono homothi a ti gips^ f^-^ 7 ^^-^^ 



reconstruction d'une image par une decomposition en 
ondelette (US - S ' 953 388 - A) , cjui eependant est 
appliqu&e la pour reconstruire una portion seulemeat de 
1» image, en exploitant la propriety de "localite" de la 
decomposition, qui n'est guere sensible aux aufcres 
portions de l' image, qu'an peut ainsi uegliger dans les 
calculs. 

L 1 invention- sera maintenant decrite plus 
ooacretetnent et conrplesfcement en liaison aux figures 
suivantea ; 

- la • figure 1 illustre une decomposition 
d r image en ondelette, 
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comprenant un serie de decompositions successives des 
jeux de projections en ondelette dormant des imagettes 
des jeux de projection, comprenant une approximation et 
des series successives de details homologues de chaque 
jeu, puis de retroprojections des imagettes des jeux de 
projection en imagettes de 1 1 image a reconstruire et 
une combinaison des imagettes a reconstruire par une 
inversion de decomposition en ondelette dormant ladite 
image, caracterisS en ce qu'il comprend, dans chacune 
des series et successivement pour des imagettes de 
details ayant un nombre croissant de points, une 
recherche de portions insignif iantes, jugees pour 
depourvues de contenu, et une recherche de portions 
homoth^tiques des portions insignif iantes pour chacune 
des imagettes de details qui suivent dans la s£rie ; et 
en ce qu ! il comprend une omission des retroprojections 
pour toutes les portions insignif iantes et toutes, les 
portions homoth^tiques . 

L 1 art anterieur comprend un exemple de 
reconstruction d'une image par une decomposition en 
ondelette (US - 5 953 388 - A) , qui cependant est 
appliquee la pour reconstruire une portion seulement de 
l f image, en exploitant la propri£t<§ de "locality" de la 
decomposition, qui n'est gu£re sensible aux autres 
portions de l 1 image, qu'on peut ainsi negliger dans les 
calculs. 

L 1 invention sera maintenant decrite plus 
concretement et completement en liaison aux figures 
suivantes : 

- la figure 1 illustre une decomposition 
d 1 image en ondelette, 
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- les figures 2/ 3 et 4 illustrent certains 
aspects de 1 1 invention, 

- la figure s illustre un organigramme 
rSsuroant true realisation complete du procede, 

5 - la figure 6 illustre un dispositif de 

prinoipe du proc£d£, 

- et les figures 7 et 8 il lust rent plus en 
detail le proc£de et cerfcaines rotations employees, 

Nous Gommexicerons par ejrposer la 

10 transformation. d l un signal par une decomposition en 
ondelette. Plusieurs modules d'ondelettes existent, qui 
ont en commun d'etre coirgparables & un filtre passe-bas* 
Le signal est s£paar£ en deux portions dont 1« une, 
associSe aux basses frequences, peut §fcre tenue pour 

15 une approximation- du signal, alors que I 1 autre , 
associ£e aux ttautes frSquenc^s, expriine plutfit aes 
details, ttoe propri^tS des ondelettes est que les 
portions peuvewt contenir chacun tine moiti£ des points 
du signal si 'bien qu»il n'y' a pas de perte 

20 deformation par cette decomposition. La decomposition 
peut etre fai fee dans le domaine direct d T expressian da 
signal ou dans le* domaine de Fourier. 

Dans ' le cas de projections d«un objefc 
examine sur un ^Sseaii g€neralement bidimensiottnel de 

25 dStecteurs, les projections peuvent Stye groupies en 
jeux bidimensioxmels d'aprSs deux de leur aoordonoSes 
(en g&n&ral p et q sur les axes d«un r£seau de 
d£tecteurs) . Toutefois, comme la figure 8 l 1 illustre # 
on considere le plus souvent des jeux; de projections 

30 r£arrang&es sur un r£seau fiotif de d6tecteurs. Dans 
1 temple de la figure 8, on regroupe des projections 
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- les figures 2, 3 et 4 illustrent certains 
aspects de I 1 invention, 

- la figure 5 illustre un organ igramme 
rlsumant une realisation complete du proc^de, 

- la figure 6 illustre un dispositif de 
principe du procede, 

- et les figures 7 et 8 illustrent plus en 
detail le procede et certaines rotations employees. 

Nous commencerons par exposer la 
transformation d'un signal par une decomposition en 
ondelette. Plusieurs modules d'ondelettes existent, qui 
ont en commun d'Stre comparables a un filtre passe-bas. 
Le signal est separ6 en deux portions dont I'une, 
associee aux basses frequences, peut etre tenue pour 
une approximation du signal, alors que l r autre, 
associee aux hautes frequences, exprime plutot ses 
details. Une propriety des ondelettes est que les 
portions peuvent contenir chacun une moitie des points 
du signal si bien qu'il n ! y a pas de perte 
d 1 information par cette decomposition. La decomposition 
peut §tre faite dans le domaine direct d 1 expression du 
signal ou dans le domaine de Fourier. 

Dans le cas de projections d'un objet 
examine sur un reseau generalement bidimensionnel de 
detecteurs, les projections peuvent etre groupees en 
jeux bidimensionnels d'apres deux de leur coordonn6es 
(en general p et q sur les axes d'un reseau de 
detecteurs). Toutefois, comme la figure 8 I 1 illustre, 
on considere le plus souvent des jeux de projections 
rearrangies sur un reseau fictif de detecteurs. Dans 
1'exemple de la figure 8, on regroupe des projections 
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ftx originaixes d f un certain nomfcre de positions Fx 
successives de la source F du rayonnemeiit, de maniere 
que les projections aboutissant & une meme colonne (S q 
constant) de d^tecteurs 5' du r§seau fictifi 4' soient 
originaires d'une taenia position Fx, et ansa! que les 
plans de projections soient tous paralldles jusqu 1 aux 
positions Ex r on a alors transforms le probl&me d'une 
georoStrie conique du, rayonneraent & une gSomStrie en 
Sventails parallSles dont la resolution et plus simple - 
De plus, le r&seau fictif 4 1 passe par le centre 0 de 
rotation et appartient done an plan de detection Pdet, 
ce qui faoilite auesi les calculs, 

L» invention pour^ait encore Stre appliquSe 
a dee rencontructions de coupes de l'objet 1 a I 1 aide 
d r u n gfegean a^nodii nen BAoni iel da d^tectertrg (toua a^ 1© 
meme coordoniiee p4} . On pourrait enc we appli^erl^ 
principe e*pos€ ci-dessus du rearrangement des 
projecl^g^jaJ^e geom&trie <kr--g?a«8<®wffla^ en fiventail 
-plaa as^gfioinatrie .parallSle. ' 

Le fcraitement des mesures r§arrang&es se 
fait etl suivant lee lignes de dgteefceure S^dvx ygae^u 
fictif 4«, successivemanc pour tous les points- 

tin signal peut Sere decompose 
succeseiveroent en ' ondelette pcur donner plusieurs 
niveau* de rSsultats, Les nouvelles decompositions ne 
concement que la portion de 1' ondelette qui donnait 
1 • approximation du signal, la portion ou lee portions 
qui donnaient les details Stant preserv&es. 

Prenons comme exemple d'objet d'une 
decomposition en ondelette une image formSe de cinq 
cercles dont un cercle ext£rieur et quatre cercles do 
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Rx originaires d'un certain nombre de positions Fx 
successivea de la source F du rayonnement, de maniere 
que les projections about issant a une meme colonne (a q 
constant) de detecteurs 5' du reseau fictif 4' soient 
originaires d'une mSme position Fx, et aussi que les 
plans de projections soient tous paralleles jusqu'aux 
positions Fx : on a alors transform^ le probleme d'une 
gSometrie conique du rayonnement a une g€om6trie en 
Sventails paralleles dont la resolution et plus simple. 
De plus, le reseau fictif 4' passe par le centre 0 de 
rotation et appartient done au plan de detection Pdet, 
ce qui facilite aussi les calculs. 

L' invention pourrait encore Stre appliqu^e 
a des rencontruct ions de coupes de l'objet 1 a 1 ' aide 
d'un reseau monodimensionnel de detecteurs (tous a la - 
raeme coordonnee p4) . On pourrait encore appliquer le 
principe expose ci-dessus du rearrangement des 
projections d'une geomdtrie en eventail plan a. 
ggometrie parallele. 

Le traitement des mesures r€arrang6es se 
fait en suivant les lignes de detecteurs 5 du reseau 
fictif 4', successivement pour tous les points. 

Un signal peut etre decompose 

successivement en ondelette pour donner plusieurs 
niveaux de resultats. Les nouvelles decompositions ne 
concernent que la portion de 1' ondelette qui donnait 
1' approximation du signal, la portion ou les portions 
qui donnaient les details etant preservees. 

Prenons comme exemple d'objet d'une 
decomposition en ondelette une image formee de cinq 
cercles dont un cercle exterieur et quatre cercles de 
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diataStres diff^reaats , toue inscribe <iaixs le premier. 
L' image da depart coznprenait n x n points, et si on 

r 

appli<jue deux fois, \ine decomposition en ondelette de 
cette image salon le prinoipe precedent, le resultafc 
S est donnS h la figure 1 „ 

ka dScomposition en ondelettes de 1* image 
dorme -un ensemble d'imagettee, dont trois sont plus 
importantes que les autres, coioprenant n/2 X n/2 points 
chacune, et correspondent aixsc details horizontaux, aiwc 

10 dStails diagonals: et aux details verticauac de 1' image 
de depart a granda Schelle ? elles sent notees 
respectivemenfc Dhl, Ddl at Dvl. details laorisontanx 

de l> image sont; obfcenus & partir des projections 
d 1 angles 6 proohe de o ou %, les detail© verticaux a 

15 partir des projection© d' angles 6 proche da jc/ 2 ou 3w/a/ 
et les details diagonaux a partir des projections 
d' angles intermediates avec les conventions de la 
figure 1. Le reste de 1' image est compose de quatre 
imagettes comprenant n/4 x n/4 point© chacune et dont 

20 trois sont das imagettes de details hori2cntaux. r 
diagonauac et verticaux 3. plus petit© echelle, notees 
©n2, Dd2 et Cv2, alors que la deroiere imagette est una 
approximation notee de 1' image de depart. Si on 
appelle ^ et y les fonctions de decomposition d'une 

25 image ou d'une imagette en ondelette, la fonction § 
donnant 1 ' approximation et la fonction y les details , 
les fonctions a appliquer a 1' image de depart pour 
obtenir la decomposition de la figure 1 sont doim€es 
par le tableau I. ' 
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diametres differents, tous inscrits dans le premier. 
L ' image de depart comprenait n x n points, et si on 
applique deux fois une decomposition en ondelette de 
cette image selon le principe precedent, le resultat 
5 est donne a la figure 1. 

La decomposition en ondelettes de 1' image 
donne un ensemble d'imagettes, dont trois sont plus 
importantes que les autres, comprenant n/2 x n/2 points 
chacune, et correspondent aux details horizontaux, aux 
10 details diagonaux et aux details verticaux de 1' image 
de depart a grande echelle ; elles sont notees 
respect ivement Dhl, Ddl et Dvl. Les details horizontaux 
de 1' image sont obtenus a partir des projections 
d' angles 6 proche de o ou %, les details verticaux a 
15 partir des projections d' angles 9 proche de it/ 2 ou 3%/ 2 , 
et les d6tails diagonaux a partir des projections 
d» angles intermediaires avec les conventions de la. 
figure 1. Le reste de 1' image est compose de quatre 
imagettes comprenant n/4 x n/4 points chacune et dont 
20 trois sont des imagettes de details horizontaux, 
diagonaux et verticaux a plus petite echelle, notees 
Dh2, Dd2 et Dv2, alors que la derniere imagette est une 
approximation no tee AA de 1' image de depart. Si on 
appelle $ et f les fonctions de decomposition d'une 
25 image ou d'une imagette en ondelette, la fonction 

dormant 1 ' approximation et la fonction \|f les details, 
les fonctions a appliquer a 1' image de depart pour 
obtenir la decomposition de la figure 1 sont donnees 
par le tableau I . 
30 AA : <t» (xi) <(> (x 2 ) <l> (^0 ♦ (i f" ) 
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Dd2 , <j> (x x ) $ (5^) ¥ (|.) ¥ (-^>, 

DVl a <j> (Xi) ¥ {Xii 
Ddl : ¥ (Xi) ¥ (Xa> 
Chi i V (Xi) + (3^) 

L' invention consists a effectues; la 
retroproj ection sur chacune dee imagettes des jeusc de 
projections decomposees en. ondelette^t a combiner les 
iro^^^^^troprojetees par une inversion da la 
jaj^^^y^ ^. en ondelette^; pour obtenir 1* image 
reeherohee, I<a decomposition en ondelettes est propice 
a diverse© simplif icatioas qui accelerent beaucoup les 
calculs, Ces siroplif icatipns sont faitee entre la 
decomposition et la cowibinaison . 

La premiere d' entre ellee concerae des 
filiations qu'on peut Stablir entre des details 
homologues a des £ciiellee diffSrentes, On considere 
pour cela des series d' imagettes dormant dee details de 
mSiue nature. La- • figure 2 (qui illustre une 
decomposition d'une • image ressemblant a celle de la 
figure 1, ma.is a •*" trois siiveaux de decomposition) 
illustre, pour les trois image ttes dee details 
horisontaux Dhl, Dhfi et Dh3, des portions homologues 
Hit «J2 et CT3 qui occupent la meme position, et la meme 
superficie relative sur chacune de ces imagettes et se 
deduisent ainsi l'une de 1' autre par une homothetie 
gSomStrique dans leurs • imagettes . 
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Dv2 : * (Xx) f (x 2 ) (J) (y) T (^>' 
Dd2 , (|> (3d) ((> (x a > T (f ) Y (^), 
Dh2 : ({) (X X ) $ (x 2 ) ¥ <^> <j) (^) 

DV1 : (J) (Xi) ¥ (X 2 ) 
Ddl : ¥ (X X ) ¥ (X 2 ) 
Dhl : Y (xx) (|) (x 2 ) 

L 1 invention consiste a effectuer la 
retroprojection sur chacune des imagettes des jeux de 
projections d£compos£es en ondelette et a combiner les 
imagettes retroprojetees par une inversion de la 
disposition en ondelette pour obtenir 1' image 
recherchee. La decomposition en ondelettes est propice 
a diverses simplifications qui accSlerent beaucoup les 
calculs. Ces simplifications sont faites entre la 
d§composition et la corabinaison. 

La premiere d' entre elles concerne des 
filiations qu'on* peut Stablir entre des details 
homologues a des gchelles differentes. On considere 
pour cela des series d 1 imagettes dormant des details de 
mime nature. La figure 2 (qui illustre une 
decomposition d'une image ressemblant a celle de la 
figure 1, mais a trois niveaux de decomposition) 
illustre, pour les trois imagettes des details 
horizontaux Dhl, Dh2 et Dh3, des portions homologues 
Jl, J2 et J3 qui occupent la meme position et la meme 
superficie relative sur chacune de ces imagettes et se 
d^duisent ainsi l'une de 1 1 autre par une homoth£tie 
geometrique dans leurs imagettes. 
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On peut fairs 1 • hypothese que, pour la 
plupart des images rencontrees en pratique (a 
1' exception notamment des images texturees) , si une 
portion telle que J3 a un contenu insignif iant, 
S c'est-a-dire qui ne revele rien par rapport atuc valeurs 
aignif icatives de 1'iroagette, les portions aontologues i. 
plus grande edhalle, ici J2 et Ji, seront elles aussi 
insignif iantes - 

Conformement a 1' invention, on. coamsnce 

10 done, pour les imagettes des d6tails, a retroproj[eter 
les details , a l'ecnelle la plus petite, puis les 
details aux echelles de plus en plus grandes. Un seuil 
numerique est applique aux coefficients de l'ondelette, 
C est -S -dire aux valeurs que prend la trans forme e dans ^ 

15 l'imagette considered „• One valeur infSrieure a ce seuil 

donne ime portion insignifiante, telle que J3» : & 
Cependant, les portions insignifiantes des imagettes ne ^ 
sent pas reconstruites, Cest>a-d±re que les calcul.s de - 
retroprojection ne sont pas effectues pour elles. 

20 On precede concretetnent <sn eonstruieant vn 

masque numerique avant de retropro j eter l»imagette. Sur 
les details horizontal Da, Is masque est coijstruif 
pour la premiere fois pour 1'imagette E>h3. II prend une 
valeur egale a 0 pour les portions insignifiantes 

25 telles que J3 et egale a 1 ailleurs. Lea coefficients 
du masque se suiveht dans un ordre determine, par 
exemple ligne apres ligne. Les calculs de 
rStroproj action sont appliques sur 1'imagette 
consideree dans 1' ordre des coefficients du masque. 

30 Quand un coefficient est egal a 0, aucun caleul n'est 
effectue pour le point correspondent de l'imagette, 
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On peut faire 1'hypothese que, pour la 
plupart des images rencontrees en pratique (& 
1' exception notamment des images texturees) , si une 
portion telle que J3 a un contenu insignif iant, 
c'est-a-dire qui ne revele rien par rapport aux valeurs 
signif icatives de 1'imagette, les portions homologues a 
plus grande 6chelle, ici J2 et Jl, seront elles aussi 
insignif iantes . 

Conformement a 1' invention, on commence 
done, pour les imagettes des details, a retroprojeter 
les details 3. l'echelle la plus petite, puis les 
details aux echelles de plus en plus grande s . Un seuil 
numSrique est applique aux coefficients de l'ondelette, 
c'est-a-dire aux valeurs que prend la transformee dans 
1'imagette consideree. Une valeur inf£rieure a ce seuil 
donne une portion insignif i ante, telle que J3 . 
Cependant, les portions insignif iantes des imagettes ne 
sont pas reconstruites , c'est-&-dire que les calculs de 
retroprojection ne sont pas effectu<§s pour elles. 

On proc^de concrdtement en construisant un 
masque numerique avant de retroprojeter 1'imagette. Sur 
les details horizontaux Dh, le masque est construit 
pour la premiere fois pour 1'imagette Dh3 . II prend une 
valeur egale & 0 pour les portions insignif iantes 
telles que J3 et ggale a 1 ailleurs. Les coefficients 
du masque se suivent dans un ordre determine, par 
exemple ligne apres ligne. Les calculs de 
retroprojection sont appliques sur 1'imagette 
consideree dans 1' ordre des coefficients du masque. 
Quand un coefficient est egal & 0, aucun calcul n'est 
effectue pour le point correspondant de 1'imagette, 
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auquel an attribue une valour nulle <?anp 1 f rmagette 
retroproject£e , 

En abordant 1' image suivance des details 
Jiori2ontaux (Dh2}, le masque numfirique permet de ne pas 
S considerer la portion J2, dont les points ne son!; pas 
trait£s par 1' unit€ de oalcvil qui efi fiectue la 
rStroprojectign. Si par example le masque do I'imagette 
Dh.3 possSde un coefficient nul 3. la ligne i et la 
coloxme j r on pr£voit que les quafcxe points des lignea 
10 2?. et 2i+i r et des colonnes 2j et 2j+l de l'iraagette 
Eh2 auronfc au©ei des valors insignif iantes- Les 
calculs d£ rStroprojection ne jseront pas effectugs pour 
ees points, 

Xto masque nmaSrique est ainsi coMtruit 

15 pour nhaque niveau de dfieornpogitioiau Le masque 
decrivant 1'itnagette Dh2 comprendra par convention de© 
coefficients egaysc S, 0 pour fcoute portion fcomologue £L 
une portion aux coefficients nuls (comma «J3) du tnasque 
de l'imagette correepondante a plus petite 6ctxelle j 

20 pour determiner les valeurs 1 ou 0 des autres points du 
masque de Db2, on utilisera encore les coinparaisoxka 
avec le seuil convent ioijnel. D'autres portions 
insignifiantes, ayant- des valeurs nulles de masque, 
ponrroAt apparaStre, C'est ce qu'on ^ illusfcre sur les 

25 . itnagettes Dvl, Dv2 * et Dv3, des decompositions 
successive** des details verticaux, L'imagette Dv3 
comprend uns portion insignifiente K3 .dont on retrouve 
les ixomologues K2 et Kl sur les decompositions a plus 
gxande ^chelle^ k'imagette Dv3 ne comprend, dans cet 

30 exemple, aucune autre portion inaignif iante, atais il a 
St«S possible d'en trouver trois autres, notSes L2, M2 
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auquel on attribue une valeur nulle dans 1 ' imagette 
retropro j ect6e . 

En abordant 1 7 image suivante des details 
horizontaux (Dh2) , le masque numerique permet de ne pas 
considerer la portion J2, dont les points ne sont pas 
traites par 1' unite de calcul qui - effectue la 
retropro ject ion* Si par exemple le masque de 1 1 imagette 
Dh3 possdde un coefficient nul a la ligne i et la 
colonne j , on pr£voit que les quatre points des lignes 
2i et 21+1, et des colonnes 2j et 2j+l de 1' imagette 
Dh2 auront aussi des valeurs insignif i antes . Les 
calculs de retroprojection ne seront pas effectu£s pour 
ces points. 

Un masque numerique est ainsi construit 
pour chaque niveau de decomposition. Le masque 
decrivant 1' imagette Dh2 comprendra par convention des 
coefficients egaux k 0 pour toute portion homologue a 
une portion aux coefficients nuls (comme J3) du masque 
de l 1 imagette correspondante & plus petite echelle ; 
pour determiner les valeurs 1 ou 0 des autres points du 
masque de Dh2, on utilisera encore les comparaisons 
avec le seuil convent i onne 1 . D' autres portions 
insignif iantes, ayant des valeurs nulles de masque, 
pourront apparaitre. C'est ce qu'on a illustre sur les 
imagettes Dvl, Dv2 et Dv3, des decompositions 
successives des details verticaux. L' imagette Dv3 
comprend une portion insignif iante K3 dont on retrouve 
les homologues K2 et Kl sur les decompositions a plus 
grande echelle. L' imagette Dv3 ne comprend, dans cet 
exemple, aucune autre portion insignif iante, mais il a 
£te possible d'en trouver trois autres, not£es L2, M2 
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et N2 f $ur l'imagette suivante JDv2„ Au niveau de 
decomposition a plus grande Schelle, celui de 
1 J imagatte Dvl, on ne s'occupes?a pas de retroprojecter 
les points situfis awe portions LI, Ml et Nl homologues 
5 de L2, M2 et K2. 

On va maintenant examiner d'autres 
particularity du procede de 1 1 invention favorables, a 
une acceleration des calculi. 

La premiere repose sur l'£galitS entre la 

10 transformSe de Fourier d'une projection de I 1 image S un 
angle fixe et la transforraSe de Fourier de 1 ' image s^r 
une droite de m§me angle passant par I'origine- 

En se reportant a la figure 3, ofi. im jeu de 
projections f'ete converti dans le domaine de Fourier A 

15 pour doxmer dee projections de nature frSquentielle 

dans le syst&&e d'aspss notSe £l et Zp r le support du i 
jeu de projection <3&ns le^d^^ine de Fourier compxrend 
las valeurs comprises entre \tffet v 0> pour k camme pour * ^ 
Une decomposition en ondelefcte, corane celle de la * 

20 figure l, fait appara£tre 1' approximation, encore notes 
RR. dans les frequences les plus basses ; autour de 
1' origin© <1 « £2 «, q # alors que les details B e 
trouvent de part, et d' autre de cette approximation./ a ' 
des frequences de plus en plus grandes pour les ditails 

25 & grande echelle. 

On a represents une droite passant par 
1' origins,, en faisanfc un angle 0 avec I'axe horizontal 
Cl. Cette droite traverse 1' approximation Azl, ainsi que 
les imagettes des details verticaux Dvl et Dv2. Mais on 
30 ' est assure que les projections faisant cet angle 0 
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et N2, sur 1'imagette suivante Dv2 . Au niveau de 
decomposition a plus grande gchelle, celui de 
1'imagette Dvl, on ne s'occupera pas de retroprojecter 
les points situes aux portions LI, Ml et Nl homologues 
5 de L2, M2 et N2 . 

On va maintenant examiner d' autres 
particularity du proced£ de 1 1 invention favorables, a 
une acceleration des calculs* 

La premiere repose sur l'egalite entre la 
10 transformee de Fourier d'une projection de 1 1 image a un 
angle fixe et la transformee de Fourier de 1' image sur 
une droite de m£me angle passant par l'origine. 

En se reportant & la figure 3, oil un jeu de 
projections a ete convert i dans le domaine de Fourier 
15 pour dormer des projections de nature f requentielle 
dans le systSme d'axes not€s C 1 et C 2 ' le support du 
jeu de projection dans le domaine de Fourier comprend 
les valeurs comprises entre -v 0 et v 0 pour £i comme pour 
£ 2 . Une decomposition en ondelette, comme celle de la 
20 figure 1, fait apparaitre 1 ' approximation encore notee 
AA dans les frequences les plus basses, autour de 
l'origine £l = C 2 = °/ alors que les details se 
trouvent de part et d' autre de cette approximation, a 
des frequences de plus en plus grandes pour les details 
25 a grande eclielle. 

On a represents une droite passant par 

l'origine, en faisant un angle 0 avec 1'axe horizontal 
£l. Cette droite traverse 1 ' approximation AA, ainsi qrie 
les imagettes des details verticaux Dvl et Dv2 . Mais on 
3 0 est assure que les projections faisant cet angle 6 




seront tout & fait inutiles pour les calculs de 
rStroprojection des details diagonawc Dd et horizontaiix 
Dh puiscjue la droite d 1 angle 0 passe a l'£cart de leurs 
ixnagettes, I*' application de 1' invention coraprend alors 
5 une selection, pour chacune des categories d'imagettes, 
des angles de projections qui seront utilises dans les 
calculs. I*e calcul sur le support des projections est 
£16mentaire «. 

La figure 4 reprend la division d'une image 

XO d£composSe en ondelette et transposSe dans le domains 
de Fourier. One reconstruction ' parf aite peut £tre 
obtenue avec 1 ' eiqploi d' un noinbre d£termin£ de 
projections, dependant de la fr£gu.ence (discretisation) 
de reconstruction de l r image de l'ofejet, Ainsi, pour 

15 reconstruire 1' approximation AA, on peut deoioatrer 
<ju'il stiff ira de choisir un noinbre de project ions s 
panpi celles qui ont £t# faites, correspondent a une 
frequence maxiTttale vl servant de rayon & un cercle 
circonscrivant la representation fr^quentielle de 

20 1< approximation &h dans la d§coroposition- De inSme, les 
details £ petite 3cbelle seront: compl&teraent rendus eo 
utilisant un norabre de projections correspondent S la 
frequence v2 dans le cercle ayan,t ce rayon , dans le 
systSme d'axes £2, et circonscrivant les 

25 representations fr£quentielles du groups de details 
Dv2, Ed2 et Dh2 a ro§me fichelle* Le support des 
projections dans le domains de Fourier permet done de 
connaitre ais&aent les frequences maximal es nScessaires 
pour la retroprojection parfaite des imagettes 

30 respective©. 
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seront tout a fait inutiles pour les calculs de 
retroprojection des details diagonaux Dd et. horizontaux 
Dh puisque la droite d 1 angle 0 passe & I'ecart de leurs 
imagettes. L' application de 1' invention comprend alors 
une selection, pour chacune des categories d' imagettes, 
des angles de projections qui seront utilises dans les 
calculs. Le calcul sur le support des projections est 
eiementaire . 

La figure 4 reprend la division d'une image 
decompos^e en ondelette et transposee dans le domaine 
de Fourier, Une reconstruction parfaite peut etre 
obtenue avec 1'emploi d'un notnbre determine de 
projections, dependant de la frequence (discretisation) 
de reconstruction de 1 1 image de l'objet. Ainsi, pour 
reconstruire 1 ' approximation AA, on peut demontrer 
qu'il suffira de choisir un nombre de projections, 
parmi celles qui ont ete faites, correspondant & une 
frequence maximale vl servant de rayon a un cercle 
circonscrivant la representation f requentielle de 
1 'approximation AA dans la decomposition. De meme, les 
details a petite echelle seront completement rendus en 
utilisant un nombre de projections correspondant a la 
frequence V2 dans le cercle ayant ce rayon, dans le 
systdme d' axes £l, Z,2 r et circonscrivant les 
representations f requentielles du groupe de details 
Dv2, Dd2 et Dh2 a m£me echelle. Le support des 
projections dans le domaine de Fourier permet done de 
connaitre aisement les frequences maximales necessaires 
pour la retroprojection parfaite des imagettes 
respectives . 
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I>& m*3e en cmvre de 1' invention pour 
exploiter ces deux particularity correprendra done, t la 
r&fcroproj action de chacune des imagettes, uae Election 
das projections utiles pour oette rStroprojection, en 
5 ecartant les autre s ; et 6venfcnellement une restriction 
du nanjbre de projections utiles ef fectiv^m^nt 
exploitSes par le calcul pour n'en garder <jue le nombre 
util^. 

I*a figure 5 est un organigrainme de 
io I'enseioble du proced£ dfecrit ici. 
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La mise en oeuvre de 1' invention pour 
exploiter ces deux particularites comprendra done, a la 
r^troprojection de chacune des imagettes, une selection 
des projections utiles pour cette retroprojection, en 
5 ecartant les autres ; et eventuelleraent une restriction 
du nombre de projections utiles ef f ectivement 
exploitees par le calcul pour n'en garder que le nombre 
utile. 

La figure 5 est un organigranutie de 
10 1' ensemble du procedg decrit ici . 
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1, Precede de reconatruction d'une image a partir de 
jeirs: de projections de cette image r comprenant une 
5 aerie de decomposition? successives des jeujc de 
projections en ondelette donnan^ele^ 
^ de3 — ^a^a-etSJtioiis , > coraprenant des images^ & f true 
approximation (AA) et de£ series successives (Dd, Dh, 
Dv) de details homologues de chaque jeu, puis dee 
10 r^troprojectiaii£s.^^^^S^3fe4&es idee jeux de projections' 0 s *? 
en imagettes de 1 1 image & reaonstruire et une 
comb inai son des imagettes & reconsfcruire par une 
inversion de decomposition dormant ladite image, 
caract6ris<§ en ce qia 5 U comprend, dans c&acune des 
IS series et success ivement pour dee imagettes de details 
ayant un norabre croissant de points , une rechezrche de 
portions insignif iantes jugSes p<yur dSpourvues de 
contenu, et une recherche de portions homoth6tiques des 1 
portions insignif iantes dans chacune des imagettes de 
20 details qui suivent dans la sfirie ; et en ce qu*il 
comprend une omission des r&troprojections pour toutes 



^^^^P^ rti ons in^ignifiatjtes efc=dSni2=^^ 

2. Precede de reconstruction d'une image gnivant^S 0 ^<hfc^ 



25 



30 



Precede de reconstruction d'une image suivant la 
revendication l, oaractSriee en ce gu»il comprend une 
selection de regions d' angles ($) des jeux de 
projections qui sont utilises dans les 
re t ropro j ect ions , 

3. Procede de reconstruction d'une image suivant la 
revendication 2, caracfc&rise en ce gu'il comprend une 
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RE VEND I CAT I ONS 



1. Precede de reconstruction d'une image a partir de 
jeux de projections de cette image,, comprenant une 
serie de decompositions successives des jeux de 
projections en ondelette dormant des imagettes des jeux 
des projections, comprenant des images d'une 
approximation (AA) et des series successives (Dd f Dh, 
Dv) de details homologues de chaque jeu, puis des 
rStroprojections des imagettes des jeux de projection 
en imagettes de 1 ! image & reconstruire et une 
combinaison des imagettes §l reconstruire par une 
inversion de decomposition dormant ladite image/ 
caracteris£ en ce qu'il comprend,' dans chacune des 
series et successivement pour des imagettes de details 
ayant un nombre croissant de points, une recherche de 
portions insignif i ant e.s jugees pour depourvues de 
contenu, et une recherche de portions homoth€tiques des 
portions insignif iantes dans chacune des imagettes de 
details qui suivent dans la serie ; et en ce qu T il 
comprend une omission des ret roproject ions pour toutes 
les portions insignif iantes et toutes les portions 
homoth£tiques . 

2. ProcSde de reconstruction d'une image suivant la 
revendication 1, caracteris<§ en ce qu'il comprend une 
selection de regions d' angles (0) des jeux de 
projections qui sont utilises dans les 
r^tropro j ect ions - 

3. Proced6 de reconstruction d'une image suivant la 
revendication 2, caract€rise en ce qu'il comprend une 
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selection d'un norabre de projections qui sont utilisees 



4. Procede de reconstruction d'une image suivant la 
revindication 2, caracterise en ce que la selection des 

5 regions d' angle est realisee en fionction d'un support 
d'une transfarmee de Fourier des jeuac de projections 

decomposes en ondeldttie. 

5, Precede de reconstruction d^une image suivant la 
revendioation 3, caracterise en ce que la selection du 

10 nombre de projections est realises en fonction de 
frequences maximales d l un support d'une transformed de 
Fourier des jeuse de proj ection^decomposes en ondelette. 
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selection d'un nombre de projections qui sont utilisees 
dans la ret roproject ions . 

4. Procede de reconstruction d'une image suivant la 
revendication 2, caracterise en ce que la selection des 

5 regions d 1 angle est r^alisee en fonction d'un support 
d'une transformee de Fourier des jeux de projections 
d£compos€s en ondelette. 

5. Procede de reconstruction d'une image suivant la 
revendication 3, caracterise en ce que la selection du 

10 notnbre de projections est r6alis£e en fonction de 
frequences maximales d'un support d'une transf orm6e de 
Fourier des jeux de projection decomposes en ondelette. 
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MESURES (JEUX BIB1MENSIONNELS 
DE PROJECTIONS, OU PREMIERE 
IMAGE DE L'OBJET) 

I 

REAMENAGEMENT DES 
MESURES CONIQUES EN 
MESURES EN EVENTAILS- 
PARALLELES OU DES 
MESURES EN EVENTAIL 
EN MESURES EN PARALLELES 



DECOMPOSITION EN ONDELETTE 
DE L' IMAGE EN IMAGETTES 



SELECTION DES 
PROJECTIONS UTILES 
POUR CHAQUE IMAGETTE 



FIG. 5 



IMAGETTES DES DETAILS A 
LA PLUS PETITE ECHELLE 



DETERMINATION D'UN 
MASQUE DES IMAGETTES 



JUL 



RETROPROJECTION 
DES IMAGETTES DE 
MEME ECHELLE 



PASSAGE A 
L'ECHELLE 
SUPERIEURE 




7V 



NON 



SYNTHESE DE L' I MAGE 
DE L'OBJET PAR COMBINAISON 
DES IMAGES RETROPROJETEES 
EN ONDELETTE INVERSE 



FIG:. 7 
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